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A precise understanding of LHC phenomenology requires the inclusion of one-loop corrections for multi-particle 
final states. In this talk we describe a semi-numerical method to compute one-loop amplitudes with many external 
particles and present first applications. 



1. Introduction 

tn this decade high energy physics will explore 
the TeV scale with hadron colliders like the Teva- 
tron and the LHC. The Tevatron experiments 
already have provided rich information on elec- 
troweak, jet and heavy quark physics. These data 
and their theoretical description are of major im- 
portance for Higgs and New Physics searches. A 
detailed understanding of Standard Model pro- 
cesses is indispensable to discriminate between 
signals and backgrounds in this respect. This is 
especially true in hadronic collisions at the TeV 
scale because of the complicated multi-particle 
final states. It is well known that leading or- 
der QCD predictions are generally plagued by 
large renormalization/factorization scale uncer- 
tainties. Only the inclusion of higher order cor- 
rections may lead to more stable predictions, as 
leading logarithmic contributions typically can- 
cel in that case. Moreover, it should be empha- 
sized that not only the overall normalization is 
an issue here, higher order corrections can also 
change the shapes of distributions. The knowl- 
edge of distribution shapes is particularly impor- 
tant in those cases where backgrounds measured 
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in a low signal region have to be extrapolated to 
signal regions by using theoretical predictions. A 
prominent example is the Higgs discovery mode 
H -> WW -> Ivl'v'. 

The computation of one-loop corrections to 
multi-particle final states is a complex task, as 
the combinatorial growth of expressions together 
with complicated denominator structures hinders 
a straightforward evaluation of the amplitudes. 
Due to its phenomenological relevance, a lot of 
activity has been going on in this direction dur- 
ing the last years |1I^I4I,5I6I7IMI0I1(>I11I1^| . 

In this talk, a method 5 for the evaluation 
of one-loop diagrams occurring in multi-particle 
computaions at NLO is presented. This approach 
allows for a good analytical control over the ex- 
pressions, but also provides stable numerical rep- 
resentations in phase space regions which typi- 
cally pose numerical problems. Some applications 
are also presented. 

2. Reduction formalism 

Any given amplitude can be represented as a 
linear combination of Feynman diagrams, which 
correspond to linear combinations of tensor inte- 
grals. We define tensor integrals by 
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where q { = k + n, r. 



and dfc = d n k/(in) n / 2 . One of the advantages 
of this representation is that the combinations 
<7, = A; + n appear naturally (e.g. fermion prop- 
agators). When introducing form factors, the 
Lorentz structure is carried by difference vectors 
A^j = rf — and g^ v . As an example, consider 
the rank 2 case 
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The kinematical information is encoded in the 
matrix 5y — (r, — Tj) — m 2 — m 2 which is the 
characteristic object in any Feynman parameter 
integral. After momentum integration, the scalar 
A-point function is of the form 
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Tensor integrals lead to additional factors 2j in 
the numerator. As a shorthand notation we use 
the ordered set S of propagator labels instead of 
the matrix S as our function arguments. A ma- 
trix corresponding to an integral with propagator 



l/Pj = 1/(<7| 
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responds to an index set S \ {j}. The reduction 
of scalar integrals can now be written as 
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One can show [S] that if ^2 bi(S) Sij = 1 is ful- 
filled, then 
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where the subscript "fin" indicates that the re- 
spective integral is infrared finite. The algorithm 
allows to isolate recursively all infrared poles in 
terms of three-point integrals. 

The reduction of tensor integrals is performed 
in a similar way. We add and subtract a linear 
combination of pinched terms and adjust the co- 
efficients conveniently: 
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proportional to a sum of higher dimensional inte- 
grals 7™ +2m (to > 0) and thus IR finite. In this 
way, the tensor integrals also can be reduced to 
a set of (potentially IR divergent) integrals with 
lower rank and less propagators and an infrared 
finite part. If S is invertible, i.e. for N < 7 and 
non-exceptional kinematics, one finds 



cl(s) = EG?- 1 ; 



otherwise the pseudo-inverse ify of the Gram ma- 
trix G can be used |5I12| to obtain 
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where the vectors Wj* span the kernel of G. 

If N > 5, one finds that all higher dimensional 
integrals drop out as a consequence of rank (S) = 
min(A, 6). Therefore one has for N > 5: 
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Note that N and the tensor rank are reduced 
at the same time. The case N — 5 is 
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more involved, but a reduction scheme has been 
worked out where no higher dimensional inte- 
grals I^ +2m (m > 0) are present for N > 4 and 
no inverse Gram determinants are introduced 5, . 
As in the scalar case, an algebraic separation of 
IR poles, contained in 3-point integrals only, is 
achieved iteratively. 

Application of the reduction formulas to inte- 
grals of rank r < 2 leads, for any N, to form 
factors in terms of simple scalar integrals. For 
rank r > 2, reduction to purely scalar integrals 
with N = 3, 4 is no longer possible without intro- 
ducing a tower of higher dimensional integrals. 
To avoid the latter, we use integrals with Feyn- 
man parameters in the numerator as reduction 
endpoints. Explicit representations for all form 
factors for r < N < 5 can be found in . The 
form factors are expressed in terms of the follow- 
ing basis integrals 
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where r max = 3 in eqs. (0) and , and of course 
purely scalar integrals 1%, J£ + 2 , i] 1 +2 , /" + 4 are 
also basis integrals. We provide two alternatives 
for the evaluation of the basis integrals. On one 
hand, an algebraic reduction to scalar integrals 
can be performed. This leads to representations 
with inverse Gram determinants. As long as the 
Gram determinant - or more precisely the related 
quantity B ~ dctG/ dct5 - is sufficiently large, 



such a representation is numerically safe. On the 
other hand, if the Gram determinants become 
smaller than a certain cut parameter A, a direct 
numerical evaluation of the building blocks of the 
reduction is performed. 

3. Numerical evaluation of basis integrals 

We show now how a numerically stable evalu- 
ation of the integrals 
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for the required cases N < 4, r < 3, d = 
4 — 2e, 6, 8 — 2e can be achieved. Let us focus 
on the basis integrals with d = 4 — 2e, N = 3 
and d = 6, N = 4. The other cases contain UV 
divergencies which are isolated in terms of exter- 
nal T(N — d/2) factors. After expansion of the 
integrand in e, the pole part is a trivial integral 
and the finite part is a logarithmic integral which 
can be treated with the same contour deformation 
methods as discussed below. IR divergencies are 
only present for N = 3 and d — 4 — 2e. For the 
IR divergent integrals explicit analytical formu- 
las can be found in The numerical evaluation 
of the remaining IR finite integrals is problem- 
atic due to kinematical singularities, which occur 
if the quadratic form x ■ S ■ x changes sign. We 
propose the following solution: We first make a 
sector decomposition 
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The integral is thus split into N terms: 

i d N (h,---,3r) = (-i) N r(N-d/2) 
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In sector I one applies the variable transformation 
Xj = tjXi for (j < I), Xj = tj-ixi for (J > I) and 
integrates out xi with the S distribution. Defining 
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Q(t) — T ■ S ■ T/2 leads to singular behaviour if 

^ N-l N-l 

Q(t) = g E + E S ^ + C = 

i,fc=i j=i 

A,B,C are denned by S. For the sector in- 
tegrals Ji the following contour deformation, 
parametrized bya,/3>0,A>0, leads to smooth 
integrands |5I14| 

x{t) = t- i 
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While AV • Q controls the size of the deformation, 
the parameters a, (3 control the smoothness of the 
deformation at the integration boundaries. There 
are situations where the contour deformation is 
not possible. This happens if V-Q — and Q = 0. 
This exceptional kinematic situation occurs in the 
presence of normal or anomalous thresholds and 
cannot be avoided. The way to deal with it is to 

split the integrations at tj = — ^2 S; if nec- 

i=i 

essary. Typically the proposed contour deforma- 
tion method works well, although the CPU time 
for evaluating a basis integral is much larger than 
evaluating the corresponding analytical represen- 
tations. The conclusion for the practitioner is 
to use the fast and accurate algebraic formulas 
for the "bulk" of the phase space and switch to 
the slow but reliable numerical evaluation near 
critical phase space regions. In the case of 4- 
point functions this looks schematically as shown 
in Fig. [U 

4. Applications 

Our reduction algorithm has been applied so 
far to several processes |15ll6ll7tT%] . Here we fo- 
cus on two of them. First, off-shell vector bo- 
son pair production via gluon fusion at the LHC, 
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Figure 1. Schematical evaluation of basic box 
functions. A is a User-defined parameter serving 
as a switch between analytic/numerical represen- 
tations. 



which is a background for the gg — > H — > W*W* 
channel |T^]. Although only box graphs occur, 
the large number of invariants (s, t, p^p^, mj, mt) 
makes the process sufficiently involved to provide 
a good testing ground for our method. The he- 
licity amplitudes were computed in a 

modular way. The amplitude was decomposed 
into gauge invariant terms and was reduced com- 
pletely to d — 6 box, d — 4 triangle and bubble 
integrals. As an illustration we show the invari- 
ant mass distribution of the charged lepton pair 
for gg -► W*W* -> Ivl'v' in Fig. H Numeri- 
cally instable regions are confined to very small 
phase space regions which do not contribute to 
the result inside the errors. 

Another application of our reduction procedure 
was the first direct computation of the gg — > jjg 
one-loop amplitude. Although the amplitude 
could be indirectly deduced from the gg — > ggg 
loop result it was shown that very com- 

pact results could be obtained for all six inde- 
pendent helicity amplitudes ^7j. This result is 
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Figure 2. Invariant mass distribution of the 
charged lepton pair for gg — > W*W* — » Ivl'v 1 . 
The two sets of curves are the LHC prediction 
with (lower) and without (upper) standard cuts. 
The effect of the third generation massive quark 
loop leads to a slight enhancement (full) com- 
pared to the case with two massless generations 
(dashed). For parameter choices see |18) . 



of relevance for the Higgs boson search channel 
PP -> H + jet -» 77 jet at the LHC. 

5. Summary 

For a precise understanding of LHC phe- 
nomenology, NLO precision for multi-particle am- 
plitudes is mandatory. We have presented a new 
semi-numerical approach for 1-loop multi-leg pro- 
cesses which is valid for an arbitrary number of 
massless as well as massive internal/external par- 
ticles. It allows to isolate IR divergences trans- 
parently and in an automated way from the am- 
plitude. Our reduction algorithm uses a set of 
building blocks which can be evaluated either an- 
alytically if the inherent Gram determinants are 
sufficiently large, or numerically, avoiding inverse 
Gram determinants and the associated instabil- 
ities completely. To achieve the latter, a multi- 
dimensional contour deformation method for 1- 
loop parameter integrals was developed. The full 
potential of the contour deformation method for 



one- and multi-loop Feynman diagrams is being 
further investigated. 

We have shown that our analytical reduction 
method works well for processes of phenomeno- 
logical interest like photon pair plus jet produc- 
tion or off-shell vector boson pair production at 
the LHC. Applications to LHC processes with a 
larger number of external particles are presently 
being studied. 
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